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ABSTRACT 

As brown dwarfs cool, a variety of species condense in their atmospheres, forming clouds. Iron and silicate 
clouds shape the emergent spectra of L dwarfs, but these clouds dissipate at the L/T transition. A variety of 
other condensates are expected to form in cooler T dwarf atmospheres. These include Cr, MnS, Na2S, ZnS, and 
KC1, but the opacity of these optically thinner clouds has not been included in previous atmosphere models. 
Here, we examine th eir effect on model T and Y dwarf atmospheres. The cloud structures and opacities are 
calculated using the Ack erman & Marlevl (1200 lb cloud model, which is coupled to an atmosphere model to 
produce atmospheric pressure-temperature profiles in radiative-convective equilibrium. We generate a suite of 
models between T e ff = 400 and 1300 K, log g=4.0 and 5.5, and condensate sedimentation efficiencies from 
/sed=2 to 5. Model spectra are compared to two red T dwarfs, Ross 45 8C and UGPS 0722-05; models that 
include clouds are found to match observed spectra significantly better than cloudless models. The emergence 
of sulfide clouds in cool atmospheres, particularly Na2S, may be a more natural explanation for the "cloudy" 
spectra of these objects, rather than the re-emergence of silicate clouds that wane at the L-to-T transition. We 
find that sulfide clouds provide a mechanism to match the near- and mid-infrared colors of observed T dwarfs. 
Our results indicate that including the opacity of condensates in T dwarf atmospheres is necessary to accurately 
determine the physical characteristics of many of the observed objects. 

Subject headings: brown dwarfs — stars: atmospheres 



1. INTRODUCTION 

Si nce the first brown dwarfs were discovered two de cades 
ago dBecklin & Zuckermanll988HNakaiima et al.1 19951) . hun- 
dreds more brown dwarfs have been discovered using wide 
field infrared surveys. These substellar objects, too low in 
mass to fuse hydrogen in their cores, range in mass from ~ 13 
to 75 Ms and are cla ssified by their spectra into L, T, and mos t 
recently Y dwarfs (iKirkpatrickl 120051; ICushing et all 120111) . 
Without hydrogen fusion as an internal energy source, brown 
dwarfs never reach a main-sequence state of constant lumi- 
nosity; instead, they cool overtime and will transition through 
the brown dwarf spectral sequence as different molecules and 
condensates form in their atmospheres. To model their at- 
mospheres accurately requires an understanding of both the 
chemistry and physics of the materials that will condense into 
clouds. 

1.1. Modeling L and T Dwarfs 
1.1.1. L dwarfs 

Grain or condensate formation has been predicted to play an 
important role in L dwarf atmosp heres since before t he firs t 
brown dwarfs were discovered dLunine et al.l 119861 [1989). 
Modern equilibrium thermochemical models predict that a 
variety of diffe rent condensates will form in L d warf atmo- 
spheres (Fegley & Lodd ers l ll994t lLoddersI 119991) : by com- 
paring models to observations, it is now well-established 
that a vari ety of refractory materials condense in L dwarfs 
(see, e.g. iTsuii et all [1991 lAllard et al.l l200!t iMarlev et al.1 
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l200l iBurrows et all 120061 ICushing et all 120081 The con- 
densates that appear to dominate, based on the abundances 
of elements available to condense, are corundum (AI2O3), 
iron (Fe), enstatite (MgSiOa), and forsterite (Mg 2 Si04), 
and these species form cloud layers, removing atoms found 
within the clouds from the lower pressure atmosph e re above 
dFeglev & Loddersl fl996|: iLodders & Feglevl [200l lLoddersI 
120031: lLodders & Fegleyll2006l: IVisscher et al.ll2010b . Within 
windows between major molecular absorption bands, there is 
little gas opacity so, in models without clouds, the emergent 
flux comes from hotter layers deep within the atmosphere. 
Cloud opacity tends to suppress the flux in the near-infrared 
within these windows; a thick cloud layer limits the depth 
from which the planet can radiate, removing some of the 
flux at these wavelengths, a nd forcing it to other wavelengths 
(I Ackerman & Marie vlBOOlT) . 

When the opacity of these clouds is included in radiative- 
convective equilibrium models of brown dwarf atmospheres, 
the res ulting model spectra match those of observed L 
dwarfs (ICushing et al.l l2006. 2008; ISaumon & Marieyll200! 
IStephens et ail2 009). Observations show that there is a range 
of colors for a given spectral type, which are believed to be 
associated with cloud variations or metallicity, but the details 
of this are not fully understood. Regardless, observed col- 
ors and spectra of L dw arfs cannot be well-m atched without a 
significant cloud layer (IBurrows et a l. 2006). 

1.1.2. T Dwarfs 

As a brown dwarf continues to cool, it undergoes a signifi- 
cant transformation in its observed spectrum when it reaches 
an effective temperature of approximately 1400 K. Objects 
cooler than this transitional effective temperature begin to 
show methane absorption features in their near-infrared spec- 
tra and, when these fea t ures appear, are cl assified as T dwarfs 
dBurgasser et al.1 120021: IKirkpatrickl 120051) . Within a small 
range of effective temperature, the iron and silicate clou ds be- 
come dramatically less important. Marl ey et al.l (12010) show 
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that this transition could potentially be explained by the break- 
ing up of these cloud layers into patchy clouds, but the de- 
tails of the transition are still very much unknown. However, 
the recent discovery of highly photometrically variable early 
T dwarfs suggests that cloud patchiness ma y indeed play a 
role (Radigan et al. 2012; Artigau et al. 2009). Regardless, as 
the clouds dissipate, the atmospheric windows in the near- 
infrared clear. Flux emerges from deeper, hotter atmospheric 
layers, and the brown dwarf becomes much bluer in J - K 
color (see Figure[TJ. 

1.1.3. History of Modeling T Dwarfs 

The first T dwarf disco v ered, G1229B (Nakaiim a et al.l 

11995^ 



19951; lOpp enheime r et al.l 



w as modeled by 
Marle v et al. l (11996b. lAllard et al.l (11996b . iFeglev & Loddersl 
(119961) . and iTsuii et alj (1 19961) using cloud-free models. 
These models assume that the condensate-forming materials 
have been removed from the gas phase, but do not contribute 
to the cloud opacity. Early T dwarfs are generally quite 
well-modeled using cloudless atmospheric models. However, 
recent observations of cooler T dwarfs suggest that T dwarfs 
of type T8 or later (T e ff< 800 K) appear to be systematically 
redder in J - K and J - H colors than the cloudless models 
predict (see Figure[TJ. 

One of the challenges of modeling brown dwarf spectra is 
the uncertainties in the absorption bands of major gas species 
such as methane and ammonia, as well as absorption due to 
collisional processes. Recent work by Sau mon et al.1 (120121) 
has modeled a rang e of brown dwarfs using improved line lists 
for ammonia from lYurchenko et al.l (1201 11) and an improved 
treatm ent of the p r essure -induced opacity of H2 collisions 
from lRichard et all (120121) . This work improves the accuracy 
of model near-infrared spectra and reddens the J-K colors 
of the model spectra with effective temperatures between 500 
and 1500 K. The color shift is due to decreased opacity in K 
band from collision-induced absorption and, for r e ff< 500 K 
model only, increased ammonia opacity in / band. However, 
these improvements do not change the colors enough to match 
the colors of the coolest T dwarfs. 

Clouds are a natural way to redden near-infrared spectra. 
Cloud opacity limits the emergent flux most prominently in 
/ ba nd, so it redden s the J -K and J-H colors of the mod- 
els. [Burgasser et al.l (120101) suggest that the remnants of the 
iron and silicate clouds could redden these cool T dwarfs, but 
here we suggest instead that the formation of other conden- 
sates, which naturally arise from equilibrium chemistry cal- 
culations, may play an important role. 

1.1.4. Y Dwarfs 

The proposed spectral class Y encompasses brown dwarfs 
that have cooled below 7^^500 K; a handful of these cool 
objects have recently b een discovered (Cushing et al.l 120111: 
Kirkpatric k et al.l2012l) . At these temperatures, NH3 begins to 
play a more significant role in shaping the near-infrared spec- 
tra, and sodium and potassium wane in importance in the opti- 
cal because they condense into clouds. Appreciable amounts 
of H2O and NH3 will condense into clouds at r e ff ^350 K and 
^200 K, respectively, and will further alter Y dwarf spectra. 
As we discover and characterize more of these cold objects, 
the study of clouds will be crucial to understand their spectral 
characteristics. 
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FIG. 1 . — Color-magnitude diagrams of L and T dwarfs. Top: Observed 
brown dwarf J —H color is plotted against the absolute H magnitude for all 
known brown dwarfs with measured parallax. M dwarfs are plotted as black 
circles, L dw arfs as red circ l es, an d T dwarfs as blue circles . Observational 
data are from Dupuv & Liu (2012). Models are plotted as solid lines. Blue 
lines are cloudless models and red lines are cloudy (/ se d=2) models that in- 
clude iron, silicate, and corundum clouds. Each labeled temperature marks 
the approximate location of the model with that effective temperature. The 
surface gravity of all models is log g =5.0 (1000m/s 2 ). Bottom: Same as 
above, but J-K color is plotted against the absolute K magnitude. 



1.2. Secondary Cloud Layers 
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Silicate, iron, and corundum, which are the condensates 
that dominate the cloud opacity in our L dwarf models, are 
not the only condensates that thermochemical models pre- 
dict will form in substellar atmospheres as they cool. Other 
condensates will form at lower temperatures and add to the 
cloud opacity via the same physical processes that formed 
the iron and silicate cloud layers. In cool substellar atmo- 
spheres, Na2S (sodium sul fide) has been predicted to play a 
poten t ially significant role (Eodders 1999; Lodders & Fegley 
2006; Visscher etaill2006T) . Other species expected to con- 
dense at these lower temperatures (roughly 600 to 1400 K) 
include Cr, MnS, ZnS, and KC1. 

To our knowledge, none of these five condensates have 
been included in a brown dwarf atmosphere model before 
now. iMarlevi (UpOO) estimated column optical depths for sev- 
eral of these species and recognized that Na2S could be im- 
portant at low r e ff but did not include this species in subse- 
quent modeling because of lack o f adequate optical constant 
data. IBurrows et all (120011 120021) noted that Na 2 S and KC1 
will condense in cool T dwarfs, but also noted that the in- 
dices of refraction are difficult to find. Helling and collabo- 
rators dHelling & Wo itke 2006) also recognized that some of 
these species will form condensates in some cases but also 
did not co mpute model at mospheres that included this opac- 
ity source. iFortnevI (120051) noted that some of these species 
might be detectable in extrasolar planet transit spectra which 
probe a longer path length through the atmosphere. 

2. METHODS 
2.1. Cloud Model 

To model cloudy T dw arf atmospheres, we modify the 
lAckerm an & Marley] (1200 ll) cloud model. This model has 
successfully been used to model the effects of the iron, 
silicate, and corundum clouds on the spectra of L dwarfs 
dSaumon & Marlevll2008HStephens et al.ll2009l) . Here, we do 
not include the opacity of iron, silicate, and corundum clouds; 
based on observed trends, we assume that the opacity of these 
clouds becomes negligible for the early T dwarfs. We instead 
includ e Cr, MnS, Na 2 S, ZnS and K C1. 

The [Ackerman & Marlev (2001) approach avoids treat- 
ing the highly uncertain microphysical processes that create 
clouds in brown dwarf and planetary atmospheres. Instead, it 
aims to balance the advection and diffusion of each species' 
vapor and condensate at each layer of the atmosphere. It bal- 
ances the upward transport of vapor and condensate by turbu- 
lent mixing in the atmosphere with the downward transport of 
condensate by sedimentation. This balance is achieved using 
the equation 



-K zz —-f sed w*q c = 0, 
dz 



(D 



where K zz is the vertical eddy diffusion coefficient, q, is the 
mixing ratio of condensate and vapor, q c is the mixing ratio 
of condensate, w» is the convective velocity scale, and / se d is 
a parameter that describes the efficiency of sedimentation in 
the atmosphere. 

This calculation provides the total amount of condensate at 
each layer of the atmosphere. The distribution of particle sizes 
at each level of the atmosphere is represented by a log-normal 
distribution in which the modal particle size is calculated us- 
ing the sedimentation flux. A high sedimentation efficiency 
parameter / se( j results in vertically thinner clouds with larger 
particle sizes, whereas a lower / se( j results in more vertically 



extended clouds with smaller particles sizes. As a result, a 
higher / se d corresponds to optically thinner clouds and a lower 
/ se d co rresponds to optically thicke r clouds. 

The lAckerman & M arley (2001) cloud model code com- 
putes the available quantity of condensible gas above the 
cloud base by comparing the local gas abundance (account- 
ing for upwards transport by mixing via K zz ) to the local con- 
densate vapor pressure /? vap . In cases where the formation of 
condensates does not proceed by homogeneous condensation 
we nevertheless comp ute an equivalent vapor pressure curve 
as described in Section l2.4.2l 

2.2. Atmosphere Model 

The cloud code is coupled to a ID atmosphere model that 
calculates the pressure-temperature profile of an atmosphere 
in radiative-conv ective equilibrium. T h e atmosphere model s 
are described in iMcKav et all (119891): iMarley et alJ (119961): 
IBurrows et"aT](ll997l):lMarlev & McKavldl999l):|Marlev et alj 
(120021) ; ISaumon & Marlevl (120081) : iFortnev et al.l (120081) . This 
methodology has been successfully applied to modeling 
brown dwarfs with both cloudy and clear atmosp heres 
dMarley et al.ll 19961 120021: IBurrows et al.ll 19971: ISaumon et al l 
1 20061 l2007pLeggett et all l2007allbl: iMainzer et al.l 120071: 
Blake et al 



I2007t ICushing et all 120081: iGeballe et ail 120091: 
IStephens et al.ll2009l) . 

In the atmosphere model, the thermal radiative transfer is 
de termined using the "source function technique" presented 
in iToon et aTl 0989). Within this method, it is possible to 
include Mie scattering of particles as an opacity source in 
each layer. Our opacity da t abase for gases, described ex- 
tensively in iFreedman et al.l (120081) . includes all the impor- 
tant absorbers in the atmosphere. T his opacity database in - 
cludes two significant u pdates since Freedman e t al.l (120081) . 
which are described in Saum on et al. (120121): a n ew molecu- 
lar line list for ammonia ( Yurchenko et al MED) and an im- 
prove d treatment of the pr essure-induced opacity of H 2 colli- 
sions dRichard et al.ll2012l) . 

Both the cloud mode l and the chemical equilibrium calcula- 
tions (see Section |2~4| | are coupled with the radiative transfer 
calculations and the pressure-temperature profile of the atmo- 
sphere; this means that a converged model will have a tem- 
perature structure that is self-consistent with the clouds and 
chemistry. 

2.3. Mie Scattering by Cloud Particles 

We calculate the effect of the model cloud distribution on 
the flux using Mie scattering theory to describe the cloud 
opacity. Assuming that particles are spherical and homoge- 
neous, we calculate the scattering and absorption coefficients 
of each species for each of the particle sizes within the model. 
In order to make these scattering calculations, we need to un- 
derstand the optical properties (the real and imaginary parts 
of the index of refraction) of each material. 

The optical properties were found from a variety of diverse 
sources, summarized in Table Q] To calculate Mie scattering 
within the model atmosphere, we use a grid of optical proper- 
ties at wavelengths from 0.268 to 227 fxm. Where data were 
not available, we extrapolated the available data, following 
trends for similar known molecules. 

The molecules with the most complete published optical 
propert i es are ZnS and KC1, both of which are obtained from 
lOuerrvl (119871) . who tabulates the optical constants for 24 dif- 
ferent minerals. 
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FIG. 2. — Na2S index of refraction. The real and imaginary parts of the 
sodium sulfide index of refraction from t he two sources used are plotted as a 
function of wavelength. Montaner et al. 11979) observational data are shown 
as a blue dashed line. Khachai et al. 1 2009) calculations are shown as a pink 
dashed line. The interpolated values used for the Mie scattering calculation 
are shown as pink circles. 



Opt ical properties for Cr are published in Stashchu k et al.l 
(1984) from 0.26 to 15 /im. The optical properties from 15 
/mi to 227 /mi were linearly extrapolated from these exper- 
imental data following the trend of the optical properties of 
Fe. Various extrapolations were tested; the choice of optical 
properties beyond 15 /im does not change the results of the 
calculations in any meaningful way. 

Optical prop e rties for MnS are published in 
Huffman & Wildl (fl96l . from 0.09 to 13 fim. Optical 
properties from 15 /im to 227 /im are extrapolated, following 
the trends of the other two studied sulfide condensates ZnS 
and Na2S. 

The optical properties for Na2S, the clouds with the largest 
optical depth, are c ombined from two different sources. 
iMontaner et al.l (119791) provid es experimental data in the in- 
frared, from 25 to 198 /jm. iKhachai et alJ ( 120091) provides 
first principles calculations of the optica l properties from 0.03 
to 91 /im. In the region of overlap, the Mont aner et al.l (11979b 
laboratory values are used. The real and imaginary parts of 
the index of refraction are plotted in Figure |2] 



TABLE 1 
Sources of Optical Properties 



Species Source 



Wavelength Range 



KC1 Ouerrv(1987) 
ZnS Ouerrv (1987) 
MnS Huffman & Wild (1967) 
Cr Stashchuk et al. (1984) 

Na 2 S Montan er et al. (1979) 
Khachai et al. (2009) 



0.22-167 /im 
0.22-167 /tm 
0.09-13 /im 
0.26-15 /tm 
25-198 /tm 
0.03-91 /tm 



TABLE 2 

Abundances of condensate-forming elements 



Element 


A(El) a 


Major condensate 


Fe 


7.54 ±0.03 


Fe 


Si 


7.61 ±0.02 


MgiSi0 4 ,MgSiCh 


Mg 


7.62 ±0.02 


MgiSi0 4 , MgSi0 3 


O 


8.76 ±0.05 


Mg^Si0 4 , MgSi0 3 , A1 2 3 , H 2 


Al 


6.54 ±0.02 


ALO3, CaAlpOi 9 , 






CaAl 2 4 , Ca 2 Al 2 Si0 7 


Na 


6.37 ±0.03 


Na 2 S 


Zn 


4.70 ±0.04 


ZnS 


Mn 


5.58 ±0.03 


MnS 


S 


7.26 ±0.03 


Na 2 S, ZnS, MnS 


Cr 


5.72±0.05 


Cr 


K 


5.18±0.05 


KC1 


CI 


5.33 ±0.06 


KC1 



NOTE. — a Where A(E1) = log[n(El) /n(H)] + 12 

2.4. Chemistry Models 
2A.1. Gas Phase Chemistry 

The abundances of molecular, atomic, and ionic species 
are calculated us ing thermochemic a l equ ilibrium f ollow- 
ing th e models of iFeelev & Loddersl (11994 119961): iLoddersI 
( 1999); Lodders & Fegl evl (120021): [Loddersl (I2OO2I 120031): 
lLodders & Fegleyl (120061) : ILoddersI d2009l) We adopt solar- 
composition element al abundances fro m lLoddersI (120031) . The 
differences betwee n ILoddersI (2003) and newer abundance 
measurements (e.g. Asplund et al. 2009) are not large enough 
to significantly a lter the conden sation temperatures consid- 
ered in the paper. Lodders (2003) abundances were therefore 
selected for consistency with previous modeling efforts by our 
groups. The abundances of condensate-forming elements are 
listed in Table 2. We assume uniform heavy element abun- 
dance ratios over a range of metallicities from [Fe/H] = -0.5 to 
[Fe/H] = +0.5 in order to explore the metallicity dependence 
of the condensation temperature expressions. 

2.4.2. Cloud Condensation Chemistry 

A simplified equilibrium condensation approach is used to 
calculate saturation vapor pressures and condensation curves 
(see Figure 3) for Cr, MnS, Na2S, ZnS, and KC1 as a function 
of pressure, temperature, and metallicity, based upon the more 
com prehen sive thermochemical models of Lodders & Fegley 
(see l2T4~TT > and Viss cher et all (120061 l2010h . In each case, we 
consider condensation from the most abundant Cr-, Mn-, Na-, 
Zn-, and K-bearing gas phases at the cloud base as predicted 
by the chemical models. The relative mass of each cloud (rel- 
ative to Na2S) is listed in Table 3, assuming complete removal 
of available condensate material from the gas phase. 

Chromium metal is the most refractory of the clouds con- 
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sidered here and condenses from monatomic Cr gas via the 
reaction 

Cr=Cr(s). (2) 

where '(s)' indicates a solid phase. The condensation condi- 
tion for Cr-metal is defined by 

P*Cr > PC, (3) 

where p' Cl is the saturation vapor pressure of Cr gas in equilib- 
rium with Cr-metal and p^. is the partial pressure of Cr below 
the cloud for a solar-composition gas (p£ r = q Cl pt, where q^. 
is the mole fraction abundance of Cr and p, is the total at- 
mospheric pressure). Upon condensation, the thermodynamic 
activity of Cr-metal is unity and the equilibrium constant (K p ) 
expression for reaction (|3) can be written as 

p'c? = K- p 1 . (4) 

Substituting for the temperature-dependent value of K p , the 
saturation vapor pressure of Cr above the cloud base can be 
estimated using the expression 

log p^r ~ 7 .490 - 20592 /T, (5) 

for T in Kelvin and p in bars. Below the cloud, we assume 
that Cr gas is approximately representative of the elemental 
Cr abundance in solar composition gas (see Table [3}: 

log p* a ~ -6.052 + log p t + [Fe/H] . (6) 

The condensation temperature as a function of the total atmo- 
spheric pressure (p,) and metallicity can therefore be approx- 
imated by setting p Cr = p' Cl and rearranging to give 

10 4 /7Lo„d(Cr) w 6.576-0.4861ogp,-0.486[Fe/H]. (7) 

This expression yields a condensation tem perature near ~ 
1520 K at 1 bar and solar metallicity (cf. Lodders & Fegley 
2006), and shows that greater total pressures and/or metallic- 
ities will lead to higher condensation temperatures. Conden- 
sation of Cr-metal effectively removes gas-phase chromium 
from the atmosphere, and the abundances of Cr-bearing gases 
rapidly decrease with altitude above the cloud. 

Our modeling of sulfide condensation chemistry follows 
that of IVisscher et al.l (12006b . and the condensation reactions 
and temperature-dependent expressions presented here are 
taken from that study. The deepest sulfide cloud expected in 
brown dwarf atmospheres is MnS, which forms via the reac- 
tion 

H 2 S+Mn = MnS(s)+H 2 , (8) 

The formation of the MnS cloud is limited by the total man- 
ganese abundance, which is 2% of the sulfur abundance in a 
solar-composition gas. The condensation curve for MnS is 
thus derived by exploring the chemistry of monatomic Mn, 
which is the domin ant Mn-bearing gas near the cloud base. 
Using results from lVisscher et all (12006 ). the saturation vapor 
pressure of Mn above the cloud is given by 

!ogPMn~ 11.532-23810/r-[Fe/H], (9) 

where the metallicity dependence comes from H 2 S (the domi- 
nant S-bearing gas) remaining in the gas phase above the MnS 
cloud base. By setting = p' Mn , th e Mn S condensation 
curve is approximated by IVisscher et al. (20061): 

10 4 /7iond(MnS)«7.447-0.421ogp,-0.84[Fe/H], (10) 



giving a condensation temperature near ~ 1 340 K at 1 bar in 
a solar-metallicity gas. 

The Na 2 S cloud is the most massive of the metal sulfide 
clouds expected to form in brown dwarf atmospheres be- 
cause Na is more abundant than either Mn or Zn in a solar- 
composition gas (see Tabled. Sodium sulfide condenses via 
the net thermoc hemic al reaction 

H 2 S + 2Na = Na 2 S(s) + H 2 . (11) 

The mass of the Na 2 S cloud is limited by the elemental abun- 
dance of sodium, which is 13% of the abun dance of sulfur 
in a so lar composition gas. Using results from Vissch eF et alj 
(2006), the saturation vapor pressure of Na above the cloud 
base is given by 

log/4, « 8.550- 13889/r-0.50[Fe/H], (12) 

where the metallicity dependence results from H 2 S remain- 
ing in the gas phase above the Na 2 S cloud and from the 
stoichiometry of Na and H 2 S in the condensation reaction. 
The condensation tem perature (where = p' Na ) is given by 
IVisscher et ail d2006l) : 

10 4 /r cond (Na 2 S) w 10.045-0.721ogp,-1.08[Fe/H], (13) 

indicating condensation near ~ 1000 K at 1 bar in a solar- 
metallicity gas. 
The ZnS cloud layer forms via the reaction 

H 2 S+Zn = ZnS(s) + H 2 , (14) 

The formation of the ZnS cloud is limited by the total Zn 
abundance, which is 0.3% of the S abundance in a solar - 
composition gas. Using results from Vissc her et al.l (2006), 
the saturation vapor pressure of Zn over condensed ZnS is 
given by 

log/4,« 12.812- 15873/r- [Fe/H] (15) 

The condensat i on cur ve (where p| n = p' Zn ) is approximated by 
IVisscher et afl(l2006l) : 

10 4 /r cond (ZnS) w 12.527-0.63 log/?,- 1.26 [Fe/H], (16) 

giving a condensation temperature of ~ 800 K at 1 bar in a 
solar-metallicity gas. 

Our treatment of KC1 condensation chemistry is similar 
to that for Cr-metal and the metal sulfides. With decreas- 
ing temperatures, KC1 replaces neutral K a s the dominant 
K-bearing gas in brown dwarf atmospheres (Loddersl ll999l 
lLodders & Feglevl 120061) . The KC1 cloud layer is thus ex- 
pected to form via the net thermochemical reaction 

KCl = KCl(s), (17) 

and condenses as a solid over the range of conditions consid- 
ered here. The vapor pressure of KC1 above condensed KCl(s) 
is given by 

log/4ci« 7.611 -11382/r, (18) 

derived from the equilibrium constant expression for the con- 
densation reaction. The mass of the KC1 cloud is limited by 
the total potassium abundance, which i s 70% of the c hlorine 
abundance in a solar-composition gas (Lodders 2003). Note 
that other K-bearing species may remain relatively abundant 
near cloud conden sation temperatures, par ticu larly at higher 
pressures (e.g., see Fegley & Lodd ersl 1 9941 and ILoddersI 1 9*991 
for a more detailed discussion of chemical speciation). How- 
ever, KC1 is the dominant K-bearing gas near the cloud base 
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TABLE 3 

Abundances of Condensate-Forming Species 



condensate 


p* below cloud base" 


cloud mass* 


Cr 


p* Cx fa 8.87 x 1(T 7 p,m 


0.30 


MnS 


p^ n « 6.32 x 10~ 7 p,m 


0.36 


Na 2 S 


p* a ?a 3.97 X Vy 6 p,m 


1.00 


ZnS 


p? Ln « 8.45 X \Q-*p,m 


0.05 


KC1 


^KCl ~ X 1" 7 Pf m 


0.12 


Fe 


pj5 e « 5.78 x 10- 5 pr'« 


20.85 


Note. — 


"Where p* is the partial pressure (in bars) 



of each gas phase species x below the predicte d cloud 
base u sing solar-composition abundances from Lodders 
12003), pi is the total atmospheric pressure, and the 
metallicity factor m is defined by logm = [Fe/H]. ''Total 
condensate mass relative to the Na 2 S cloud. Values for 
Fe shown for comparison. 

for the relevant P— T conditions expected in cool brown dwarf 
atmospheres (see Figure 3) over the range of metallicities (- 
0.5 to +0.5 dex) considered here. For simplicity we therefore 
assume that KC1 is approximately representative of the ele- 
mental K abundance below the cloud, given by 

logPKCi « -6.593 + logp, + [Fe/H] . (19) 

The condensation temperature as a function of pressure and 
metallicity is estimated by setting p^ a = p' KCl and rearranging 
to give 

10 4 /7iond(KCl) » 12.479-0.8791ogp,-0.879[Fe/H], (20) 

yielding a condensation te mperature near ~ 800 K at 1 bar in 
a sola r-metallicity gas (cf. Lodders 1999; Lodde rs~& Feglevl 
2006]). In general, the condensation curve expressions demon- 
strate that condensation temperatures increase with total pres- 
sure, as illustrated in Figure 4. Furthermore, higher metallici- 
ties are expected to result in higher condensation temperatures 
and more massive cloud layers in brown dwarf atmospheres. 
In each case, the saturation vapor pressures of cloud-forming 
species rapidly decrease with altitude above the cloud layers. 

2.5. Comparison to Other Cloud Models 

The |Ackerman & Marley] (120011) model is one method of 
several t hat have been applie d to cloudy brown dwarf atmo- 
spheres. Hellin get al.l d2008) review various cloud modeling 
techniques and compare model predictions for various cases. 
The most important conceptual differences between these ap- 
proaches lies in the assumptions of how condensed phases in- 
teract with the gas. 

In the chemical equilibrium approach (e.g. Allard et al. 
120011) . condensed phases remain in contact with the gas phase 
and can continue to react with the gas even at temperatures 
well below the condensation temperature. As an example, 
when following this approach, Fe grains which first condense 
at temperatures of over 2000 K react with atmospheric H2S 
to form FeS below 1000 K. In the condensation chemistry ap- 
proach we employ here, the condensed phases are assumed 
to sediment out of the atmosphere and are not available to 
interact with gas phases at temperatures below the conden- 
sation temperature. Thus Fe grains form a discrete cloud 
layer and do not react to form FeS. H2S consequently re- 
mains in the gas phase and reacts to form condensates as 
outlined in Section 2.4.2. Jupiter is an excellent example of 
the applicability of this framework, as the presence of H2S 



in the observable atmosphere is only possible because Fe is 
sequestered in a deep cloud layer, which prevents the forma- 
tion of FeS which otherw ise deplete other gas phase S species 
dFegley & Lodderslll99l . The presence of alkali absorption 
in T dwarfs likewise demonstrates the applicability of con- 
densation chemistry (Ma rley et al.l l2002). A detailed compar- 
ison of true equilibrium condensation and cloud condensate 
remov a l from equilibrium can be found in Feglev & Lodders 
dl994l) ; lLodders & Feglevl (12006b and references therein. 

A different approach is taken by Helling & Woitke 
dHelling & Woitkel 12006) who follow the trajectory of tiny 
seed particles of TiC>2 that are assumed to be emplaced high 
in the atmosphere and sink downwards. As the seeds fall 
through the atmos p here t hey collect condensate material. In 
Helling & Woitke (2006) and numerous follo w on papers 
dHelling et all l2008t IWitte et alJ 120091 1201 U Ide Kok et alJ 
1201 ll) this group models the microphysics of grain growth 
given these conditions. Because the background atmosphere 
is not depleted of gaseous species until the grains fall through 
the atmosphere, a compositionally very different set of grains 
are formed. In particular they predict 'dirty' grains composed 
of layers of varying condensates. 

A direct comparison between the predictions of the various 
cloud modeling schools is often difficult because of differ- 
ing assumptions of elemental abundances and the background 
thermal profile. Modeling tests in which predictions of the 
various groups are compared to data would be illuminating, 
but this is far beyond the scope of the work reported here. 

2.6. Evolution Model 

In order to calculate absolute magnitudes of the modeled 
brown dwarfs, we use the results of evolution models which 
determine the radius of a brown dwarf as it cools and con- 
tracts over its life t ime. We use the evolution models of 
Saumon & Marley (2008) with the surface boundary condi- 
tion from cloudless atmospheres. Using a cloudless boundary 
condition instead of one consistent with these clouds changes 
the calculated magnitudes of the models very slightly, but 
does not change the overall trends or results. 

2.7. Model Grid 

To analyze the effect of these clouds, we generate a grid of 
182 model atmospheres at effective temperatures and surface 
gravities spanning the full range of T dwarfs. We calculate 
pressure-temperature profiles and synthetic spectra for atmo- 
spheres from 400 to 1300 K (50 to 100 K increments), with 
log(g) (cgs) of 4.0, 4.5, 5.0, and 5.5 and cloud sedimenta- 
tion efficiency parameter / se d=2, 3, 4, and 5. For this study, 
we use only solar metallicity composition. We then compare 
these, both photometrically and spectroscopically, to observed 
T dwarfs. 

3. RESULTS 

3.1. Model Pressure-Temperature Profiles 

In Figure [3] we show the pressure-temperature profiles of 
models with effective temperatures of 400 K, 600 K, 900 K, 
and 1300 K. The surface gravity of the 400 K models is log 
g=4.5; for the hotter models, log g=5.Q. We plot models with 
two different cloud sedimentation efficiencies which include 
only the Na 2 S, MnS, ZnS, Cr, and KC1 clouds. Because Na 2 S 
and MnS are by far the most dominant cloud species (see Sec- 
tion |3]3]l, as a shorthand we refer to this collection of clouds 
as 'sulfide clouds.' 
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Temperature (K) 

FIG. 3. — The pressure-temperature profiles of model atmospheres are plot- 
ted. Models at 400, 600, 900, and 1300 K are shown, and the effective tem- 
perature of the model is labeled on the plot. The surface gravity of the 400 
K model is log g=4.5; for the hotter models, log g=5.0. We show cloudless 
models in blue, and cloudy models with / se d=2 (red) and 4 (orange). Con- 
densation curves for each condensate species are plotted. The cloudy models 
include the condensates Cr, MnS, Na2S, ZnS, and KC1. Note that for each 
case, increasing the cloud thickness increases the temperature at a given at- 
mospheric pressure. The 1-6 fim photosphere of each model is shown as a 
thicker line. 

The condensation curves of major and minor species pre- 
dicted to form by equilibrium chemistry are also plotted. The 
location of a given cloud base is expected to be where the 
pressure-temperature profile of the model atmosphere crosses 
the condensation curve. Each model crosses the condensa- 
tion curve of each species at very different pressures and, to a 
lesser extent, temperatures, so we expect that the significance 
of the clouds will be strongly controlled by the effective tem- 
perature of the model. 

The cloudy 400 K and 600 K models have two convection 
zones. All 900 and 1300 K models have a single deep convec- 
tion zone. 

For all models, it is clear that, as in previous cloudy models 
of L dwarfs, adding cloud opacity has a "blanketing" effect 
on the model, increasing the temperature of the atmosphere 
for a given atmospheric pressure. As the cloud becomes opti- 
cally thicker, the entire pressure-temperature profile becomes 
hotter; thus, on a plot of pressure-temperature profiles such as 
Figure[3] increasing the cloud opacity moves the whole profile 
to the right. 

3.2. Model Spectra 

In Figures [4] and [5] we show the spectra of the same ex- 
ample models at 1300, 900, 600, and 400 K: Figure shows 
the wavelength-dependent brightness temperatures from these 
models, while [5] shows the model fluxes computed from the 
top of the atmosphere. The brightness temperature gives some 
insight into the depth into the atmosphere probed at each 
wavelength. Flux from 0.8 to 1.3 /im comes from the deep- 
est, hottest layers of the atmosphere. Clouds change the flux 
in this wavelength range by limiting the depth from which 
flux emerges; conversely, the clouds do not change the depth 
probed between ^2 and 5 /im because the clouds form below 
the layers from which most of the flux is emerging. However, 
the hotter atmospheric temperatures at a given pressure (see 
Figure 3) lead to slightly higher fluxes at these wavelengths. 
Though not plotted here, flux in the mid-infrared also comes 



from above the cloud layer and is slightly higher because the 
entire pressure-temperature profile is hotter. 

For the hottest of these models (T e ff=1300 K), the cloudy 
spectra look almost identical to the cloudless spectrum. This 
model is too hot to have much mass of these condensed 
species form in the photosphere. For a cooler model (T e ff=900 
K), the cloudy spectra look different from the cloudless spec- 
trum. As we decrease the sedimentation efficiency / se d in the 
model, the flux in Y and J bands decreases and the flux in K 
band increases. 

For the coldest two models shown (T e ff=400, 600 K), the 
cloudy spectra look dramatically different from the cloudless 
spectrum in the near-infrared; even the thinnest cloud consid- 
ered here (/ se d=5) causes the flux in Y and J to decrease by 
50% and the flux in K to correspondingly increase. Decreas- 
ing the sedimentation efficiency enhances this effect. 

Figure [6] shows the effect of the iron and silicate clouds on 
the spectra of models with the same effective temperatures 
(1300 K, 900 K, 600, and 400 K) and surface gravity. Note 
that unlike the sulfide clouds, these iron and silicate clouds 
substantially change the shape of the 1300 K and 900 K mod- 
els by suppressing the flux in Y, J, and H and increasing the 
flux in K band and the mid-infrared. This strong effect at 
higher temperatures is due to the large amount of iron and 
silicate condensed in the visible atmosphere at those temper- 
atures. 

3.3. Cloud Structure in Model Atmospheres 

Figure [7] shows the distribution of clouds in the model at- 
mospheres for the three example cases. The locations of iron, 
silicate, and corundum clou ds, using models tha t only i nclude 
those clouds — the standard Saumon & Marlevl d2008b cloud 
configuration — are plotted for r eference. The column optica l 
depth is given by Equation 16 in[Ackerman & Marlevl (2001). 
which calculates the cumulative geometric scattering optical 
depth by cloud particles through the atmosphere. 

For the 1300 K model, all of the sulfide clouds have tiny 
optical depths in the photosphere and do not affect the emer- 
gent spectra. The silicate and iron clouds would have signifi- 
cant optical depth (r=2-3) and would substantially change the 
emergent spectra. 

For the 900 K model, all of the sulfide clouds have a col- 
umn optical depth smaller than 1 in the photosphere. KC1 and 
ZnS have tiny optical depth (r < 2 x 10~ 2 ) and will not create 
an observable change in the spectrum. Na2S and MnS have 
optical depth between 0. 1 and 1 and will change the model 
spectra slightly. 

For the 600 K model, Na2S is the most important conden- 
sate opacity source. KC1 has a small optical depth and ZnS 
has a negligible optical depth. This result is expected, based 
on the abundances of each species (see Table [3). The other 
two clouds, MnS and Cr, are below the near-infrared photo- 
sphere, so also do not change the spectrum. The silicate and 
iron clouds would also be below the photosphere. 

Using our full grid of models, we can examine the impor- 
tance of the Na2S cloud as a function of T & s and surface grav- 
ity. Figure|8]shows how the column optical depth of this cloud 
varies from 400-1300 K and log g from 4.0-5.0 for a constant 
value of the parameter / sec |. Moving to r e ff values below 1300 
K, the Na2S cloud grows in importance as it forms progres- 
sively deeper in the atmosphere, so that there is larger mass 
of condensate in the cloud. The maximum optical depth of 
the Na2S cloud at pressure levels above the bottom of the 1- 
6 iim photosphere is largest for models with T e s of 600 K. 
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FIG. 4. — The model spectra are plotted as brightness temperature vs. wavelength, cloudless, / se d=2, and / se d=4 models are shown. The solid horizontal line 
indicates the temperature at the base of the each cloud, and the dashed horizontal line denotes the temperature of the layer in which column extinction optical 
depth of the cloud reaches 0.1. Note that for all clouds in the 1300 K model, the column optical depth model never exceeds 0.1. The maximum column 
optical depth of the Na2S clouds (t at the cloud base) is calculated using the / scc i=4 models and labeled on each plot. 



At lower r e ff, much of the cloud opacity is below the visible 
atmosphere. The optical depth within the photosphere is sig- 
nificant (t > 1) for models between 400 and 700 K. 

Figure [8] indicates that the cloud optical depth within the 
photosphere is largest for higher surface gravity atmospheres 
for a constant value of / se( j. However, this does not necessar- 
ily predict that higher gravity brown dwarfs will have thicker 
clouds than lower gravity brown dwarfs because the parame- 
ter / S ed is not necessarily independent of gravity. 

4. COMPARISON WITH OBSERVATIONS 

4. 1 . Color-Magnitude Diagrams 

In Figure [9] we plot the photometric colors in the near- 
infrared of all brown dwarfs with measured parallaxes and 
appar ent J and K magnitude errors smaller t han 0.2 magni- 
tudes (Du puv & LiullioialFahertv et al.ll2012l) . We also plot 
the calculated photometric colors of our suite of cloudless 
and cloudy models from 400-130 K with several represen- 
tative surface gravities. In Section |4~T1 we discuss the general 



trends of the photometric colors of models at various effective 
tem peratures and cloud sedimentation efficiencies. In Section 
14.21 we compare our model results to the photometric obser- 
vations. 

As discussed in Sections [3. 21 and [331 clouds in T dwarf at- 
mospheres tend to suppress the flux in Y and J bands and in- 
crease the flux in K band. The flux shift from J to K gives 
cloudy models larger (redder) J — K and J-H colors than 
cloudless models. 

In Figure|9] the hottest cloudy models have nearly the same 
near-infrared colors as cloudless models. As we decrease the 
effective temperature of a cloudy model, more cloud material 
condenses; the model has a redder photometric color than a 
cloudless model with the same 7^. If we reduce the sedimen- 
tation efficiency (/ se d) of the cloud, the cloud becomes opti- 
cally thicker, and the model has a redder photometric color. 

The upper panels, which show J — K photometric colors, 
show that our sulfide cloud models can easily reach the colors 
of red T dwarfs, with / se d values of 4-5. The bulk of the T 
dwarf population is bluer than the / se d=5 model. However, the 
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FIG. 5. — Model spectra. From top to bottom, T cf f=1300 K, log g=5.0; 
T c ff=900 K, log g=5.0; r eff =600 K, log g=5.0; r eff =400 K, log g=4.5. We 
show cloudy models with / sc d=2 and 4 which include the condensates Cr, 
MnS, Na2S, ZnS, and KC1 and cloudless models for comparison. Note that 
for the r e ff=400 K, T c ff=600 K and 7^=900 K models, the cloudy models are 
progressively fainter in Y and J bands and brighter in K band as the sedimen- 
tation efficiency decreases. In contrast, for the r e ff=1300 K case, the clouds 
do not significantly change the spectrum. 
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FIG. 6. — Model spectra with iron/silicate clouds. As in Figure [5] from 
top to bottom, T cff =1300 K, log g=5.0; r eff =900 K, log g=5.0; r cff =600 
K, log g=5.0; T e ff=400 K, log g=4.5. We show cloudy models with 
iron/silicate/corundum clouds (no sulfide clouds) with f se( i=2 and cloudless 
models for comparison. Note that these clouds, unlike the sulfide clouds in 
Figure|5] significantly change the shape of the 1300 K model. 
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FIG. 7. — Pressure vs. column optical depth. The column optical depth of 
each cloud species is plotted. The solid lines denote the clouds examined in 
this study: NaiS, ZnS, KCI, Cr, and MnS. The dashed lines show the column 
optical depths of models that include only the more refractory clouds corun- 
dum (AI9O3) iron (Fe), and forsterite (Mg2Si04) to show where those clouds 
would form in comparison to the sulfide clouds. All models use /sed-2- The 
shaded grey area shows the region of the atmosphere which lies within the A = 
1 to 6 fim photosphere. Note that the Na2S cloud is by far the most important 
of the added clouds for the 600 K model in the near-infrared. Also note that if 
the AI2O3, Fe, and Mg9Si04 persisted to effective temperatures of 900-1300 
K, they would be quite visible, which would not match observations. 



I 1 1 1 1 I 1 1 1 1 I 1 1 1 
•— • log g: 5.0 
•— • log g: 4.5 
•— • log g: 4.0- 




600 700 800 900 1000 1100 1200 1300 
Model effective temperature (K) 

FIG. 8. — The column optical depth of the Na2S cloud above the bottom 
of the 1-6 fim photosphere in each model is plotted as a function of model 
effective temperature. The curves connecting the points are there to guide the 
eye. Three different surface gravities are shown and all models use / S ed=2. 
The column optical depth peaks at temperatures of about 600 K, and models 
with higher surface gravity have a greater Na2 S column optical depth. 



cooler T dwarfs are generally well-matched by the models. In 
J-H, the color directly affected by cloud opacity limiting 
the depth to which one sees, the model colors are an excellent 
match to the data. The cloudy models are a much better match 
than the corresponding cloud-free models. 

4.2. Comparison to Observed T Dwarfs 
4.2.1. Expected Surface Gravity of T dwarf s 

Based on Sau mon & Marie vl (120081) evolution models and 
assuming that observed brown dwarfs will have ages less than 
10 Gyr, we expect that the coldest objects modeled, between 
400-600 K, will have surface gravities less than log g=5.2. 
Hotter objects, between 1000 and 1300 K, will have surface 
gravities less than log g=5.5. 

4.2.2. Cloud Sedimentation Efficiency 

For the L dwarfs, we are generally able to match pho- 
tometric colors by including silicate, iron, and corundum 
clouds with a sedimentation efficiency paramet er of / se d=2±l 
dStephens et al.l2009HSaumon & Marlevl2008l) . However, for 
these cooler T dwarfs, models with these sulfide clouds with 
/ se d=2 are redder than observed brown dwarfs. If we assume 
/ S ed is larger — around 4 or 5 — we are able to match observed 
colors. The cloud model does not explicitly suggest any phys- 
ical mechanism for why / se d would be different. However, 
since these objects are about 1000 K cooler than L dwarfs, 
it would not be surprising if these objects populate a differ- 
ent physical regime, and would have substantially different 
rates of atmospheric mixing and cloud condensation. Indeed, 
a large increase in / se d with lower T e s values i s one way to 
quick ly clear away the silicate and iron clouds (Kna ppet alj 
12001 

4.2.3. WISE Color-Color Diagrams 

ICushing et alj d201 lh announced the discovery of six pro- 
posed Y dwarfs found by the WISE mission. Objects around 
the T-to-Y transition are cold enough to have NH3 absorption 
features in the near infrared. 
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FIG. 9. — Color-magnitude diagrams for M, L, and T dwarfs. As in Figure[T] observed ultracool dwarf color is plotted against the absolute magnitude for all 
known brown dwarfs with measured parallax. In the top 3 plots, J-K color is plotted against absolute J magnitude; in the bottom 3 plots, J-H color is plotted 
against absolute H magnitu de. All photometry is i n the MKO system . M dwarfs are plotted as black circles, L dwarfs as red circles, and T dwarfs as blue circles. 
Observational data are from Dupuv & Liu (2012); Faherty et al. 12012). The locations of the brown dwarfs Ross 458C and UGPS 0722-05, the objects to which 
we compare model spectra to observations in Figures 11 and 12, are shown with a purple star and square symbol, respectively. Models. Models are plotted as 
lines. Each labeled temperature marks the approximate locations of the model with that effective temperature. Three representative gravities are plotted: from 
left plot to right plot, log g=4.0, 4.5, and 5.0. Blue lines are cloudless models and red lines are cloudy models (/ sc( ] = 5, 4, 3, and 2 from left line to right line in 
each plot) that include the opacity of only the newly added clouds — Na2S, Cr, MnS, ZnS, and KC1. 



We have obtained near-infrared photometry for two pro- 
posed Y0 dwarfs discovered using the Wide-field Infrared 
Survey Explorer (WISE) by Cushing et al. (2011). YJH 
was obtained for WISEP J140518.40+553421.5 and YJ for 
J15415 1.65-225025 . 2, using the near-infrared imager NIRI 
dHodapp et al.l 120031) on the Gemini-North telescope on 
Mauna Kea, Hawaii. T he photometry is on the M auna Kea 
Observatories system (Tokunaga & Vaccal 120051) . The data 
were obtained via the Gemini queue program GN-2012A-Q- 
106 on 2012 February 10; the program aims to supplement 
and improve on the photometry presented in the discovery 
paper. Individual exposure times of 60 s were used at Y and 
J, and 30 s at H; a 9-position dither pattern with 10" offsets 
was repeated as necessary for sufficient signal to noise. 
The total exposure time for WISEP J140518.40+553421.5 
was 9 minutes at Y and J and 58.5 minutes at H\ for 
WISEP J154151.65-225025.2 the total exposure time was 
18 minutes at each of Y and J. Data were reduced in a 
standard fashion and flatfielded with calibration lamps on the 
telescope. The UKIRT faint standards FS 133 and 136 were 



used for photometri c calibration; J and H were taken from 
Leggettetal. (2006), and Y from the UKIRT online catalog 
(http://www.jach.hawaii.edu/UKIRT/astronomy/calib/phot 
cal/fs ZY MKO wfcam.dat). 

The final reduced magnitudes are: Y = 21.41 ±0.10, 
J = 21.06 ± 0.06 and H = 21 Al ± 0.08 for WISEP 
J140518.40+553421.5; 7 = 21.63 ±0.13 and / = 21. 12±0.06 
for WISEP J15415 1.65-225025.2. 

Figure [10] shows how clouds will effect these cold objects. 
H — W2 is a useful temperature indicator for these objects, 
while J-H is sensitive to both the cloud structure and grav- 
ity. As the r e ff of the non cloudy models decreases from 800 to 
500 K, the models become progressively bluer in J-H color. 
However, most of the proposed WISE Y dwarfs are redder 
than the cloudless model. The models that include the sulfide 
clouds match their colors better. This result tentatively sug- 
gests that for objects colder than T dwarfs, the sulfide clouds 
remain important. Of course, for objects with effective tem- 
peratures of ^350 K, water will condense; at t hat point, H2O 
clouds should contribute to the spectra (e.g. iBurrows et al.l 
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FIG. 10. — Color-color diagra ms using WISE and near infrared data . Observed J—H versus H — W2 colors of L and T dwarfs (Kirkpatrick et al. 2011) and 
proposed WISE Y dwarfs (Cushing et al. 201 1 ; Kirkpatrick et al. 2012) are plotted. For J and H bands we use MKO photometry. L and T dwarfs are plotted as 
red and blue dots, respectively. WISE Y dwarfs are plotted as purple error bars; Y dwarfs with magnitude upper limits are shown in pink. Model photometric 
colors are shown as solid and dashed lines; the blue line shows a cloudless model and the red lines show two cloudy models (from left to right, / sc d=4 and / sc d=2). 
Each labeled temperature marks the approximate location of the models with that effective temperature. Many of these cold brown dwarfs have photometric 
colors closer to the cloudy models than the cloud-free model. The left plot shows log £=4.0 and the right plot shows log g=5.0. 



2003). 

4.3. Comparison to Observed T Dwarf Spectra 

We now compare model spectra to two relatively cold, red 
T dwarfs, Ross 458C and UGPS 0722-05. The near-infrared 
spectra of these two objects are not well-matched by cloudless 
T dwarf spectra; by including our neglected clouds, which for 
these cool objects are dominated by the Na2S cloud, we match 
their spectra more accurately. 

We compare models to both near-infrared spectra and near- 
and mid -infrared photometry . As in previous studies of brown 
dwarfs dCushing et al.1 12001 iStephens et al] [20091) . we find 
that in different bands, the observations are best fit by mod- 
els of different parameters. In this study, we focus on finding 
models that fit the shape of the spectra in the near-infrared 
where clouds play a significant role. 

4.3.1. Ross458C 

Ross 458C is a late-type T dwarf (T7-9) which is sepa- 
rated by over 1 100 AU from a pair of M star primaries. It has 
anomalously red near- infrared colors (J-K = -0.21 ± 0.06). 
iBurgasser et al.l (120101) obta ined spectroscopic obser vations 
with the FIRE spectrograph (ISimcoe et al.ll2008l |2010[) on the 
Magellan Baade 6.5 meter telescope at Las Campanas Obser- 
vatory. They fit the spectrum using cloudless and cloudy mod- 
els (which include only the opacity of the iron, silicate, and 



corundum clouds) and find that cloudy models fit significantly 
better than cloudless models. Burgasser et al. (2010j) con- 
clude that cloud opacity is necessary to reproduce the spec- 
tral data and invoke a reemergence of the iron and silicate 
clouds. We instead assume that the iron and silicate clouds 
are depleted, as we observe generally in other T dwarfs, and 
investigate the effect of the sulfide clouds. 

In Figure QT| we show the FIRE spectrum and the best 
fitting cloudy and cloudl ess models. We al so show pho- 
tometry i n J, H, and K dDupuv & Liul 120121) , WISE pho- 
tometry (iKirkpatrick et al.l 120121) . and Spitzer photometry 
dBurningham et alj|201 ll) . The spectr a used to generate these 
results differ somewhat from those in Burg asser et al.l (120101) 
because we use mod els that include recen t improvements to 
the opacity database (ISaumon et al.ll2012l) for both ammonia 
and the pressure-induced opacity of H2 collisions. All models 
are fit by eye to the observatio ns. 

Like Burgasse r et al.l (120101) . we find that clouds are essen- 
tial to match the spectrum of Ross 458C. FigureQT]shows the 
best fitting cloudless model and the two best fitting cloudy 
models (one including the iron and silicate clouds and the 
other including the sulfide clouds). The cloudless model is 
a poor representation of the spectral data; the flux in Y and / 
is too high and the flux in H and K is too low. The cloudy 
models are better representations of the relative flux in each 
band. 
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FIG. 1 1 . — Ross 458C near-infrare d spectrum comp arison between data and models. Three different models are compared to the observed spectrum and 
photometry of Ross 458C from Burgasser et al. (2010). The left panels show the near-infrared spectra; the right panels show the same spectra and models with 
near- and mid-infrared photometry. Yellow points show J, H, and K photometry; orange show Spitzer [3.6] and [4.5] photometry; red show WISE Wl, W2, 
and W3 photometry. The filters for the photometric bandpasses are shown with corresponding colors along the bottom. The top row shows a cloudless model 
spectrum that best matches the data, which has an effective temperature of 800 K and surface gravity log g=4.0. The middle row shows the best matching cloudy 
spectrum using iron and silicate clouds. The bottom row shows the best matching cloudy spectrum using sulfide clouds. Both cloudy models have significantly 
lower effective temperature (100-250 K cooler) than the cloudless best-matching model, but have similar (low) surface gravity. For this object, we can largely 
match the overall features of the spectrum using either iron/silicate clouds or sulfide clouds. 
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Burgasse r et al.l (120101) found that the surface gravity of 
Ross 458C must be low (log g=4.0) for models to match the 
observed spectrum. Likewise, we find that our best-fitting 
models have surface gravities of 4.0 (cloudless), 3.7 (silicate 
clouds), and 4.0 (sulfide clouds). 

We conclude that we do not need to invoke a reemergence 
of iron and silicate clouds into the photosphere of Ross 458C 
to reproduce the observed spectrum. Instead, we are able to 
reproduce the spectrum using the sulfide clouds which are nat- 
urally e xpec ted to form in the photospheres of cool T dwarfs. 
Section 15.21 contains additional discussion on which cloud 
species we expect to be important. 

The very red slope to the L band spectrum of Ross 458C — 
much redder than all the models — is reminiscent of the be- 
havior of some cloudy L dwarfs, including 2MASS2224 and 
DE 0255 (L3.5 and L9 r espectively) and may be a signature 
of very small dust grains (Stephens etaDl2009h . 

The discrepancies at 4.5 /im are likely to be a result of non- 
equilibrium chemistry, which is not included in t hese models. 
This effect is discussed in more detail in Section l4~4l 

4.3.2. UGPS J072227. 5 1-05403 1.2 

UGPS J072227.5 1-05403 1.2 (hereafter UGPS 0722-05) is 
a T9 or T10 dwarf with an ef fective temperature of approx- 
imately 500 K, discovered by iLucas et aD (120101) . Previous 
spectral analysis with cloudless models has been unsuccess- 
ful at modeling the flux in the near-infrared in Y and J bands 
dLeggettet al.ll2012l) . 

In F igure fT2l we plot the near-infrared spe ctra published in 
ILucas et all (120101) and lLeggettet al.l (120 121) with the cloudy 
and cloudless models that are fit by eye to be the closest 
representations of t he data. We als o show J, H, K and 
Spitzer photometry (ILucas et al.ll20 1 Qb and WISE photometry 
(Kirkpatrick et al. 2012|). These models have similar temper- 
atures and gravities to previous studies; Legg ettet al.l (120121) 
presented fits with T e s between 492 and 550 K and log g=3.52 
to 5.0, whereas our fits have T e ff of 600 K (cloudless) and 500 
K (with sulfide clouds) and both have log g=4.5. 

Note that the flux in Y and J in cloudless models is system- 
atically too high. The opacity of the sulfide clouds provides 
a natural mechanism to decrease the flux in Y and J and in- 
crease the flux in H and K bands. The match to the models 
is still not perfect. This may be due in part to incomplete line 
lists for methane; these cold objects have a significant amount 
of CH4 in the atmosphere, which absorbs strongly in the near 
infrared. The discrepancy at 4.5 fim is most likely due to ab- 
sorptio n of CO as a result of non-equilibrium chemistry (see 
Section l4~4l i. Out stand ing issues in T dwarf modeling are dis- 
cussed in Section IB31 

4.4. Non-Equilibrium Chemistry 

We note that both of our preferred sulfide cloud models 
predict brighter M band (4.5 fim) photometry than is ob- 
served. This is likely a consequence of our neglect of non- 
equilibrium mixing of CO in this study. As first predicted by 
Fegley & Lodders (1996), such mixing is a n important pro- 
cess in the atmosp heres of brown dw arfs (|Noll et al.l 119971: 
iMarley et aD2002l: ISaumon et al.ll2006l IStephens et al.ll2009l) 
as it is in solar system giants (Bars hay & Lewisll 19781) and the 
rela tive impact of the mixing increases with decre asing grav- 
ity dHubenv & Burrowsll2"007t iBarman et al]|201 ll) . Absorp- 
tion by excess atmospheric CO decreases the thermal emis- 
sion in M band and is likely responsible for the mismatches 



seen in Figures QT| and Q~2] particularly for the lower gravity 
Ross458C. 

The formation of the clouds considered in this study does 
not involve the species most affected by non-equilibrium 
chemistry such as CO and CH4. The cloud models will 
therefore be only minimally affected by the changes in the 
pressure-temperature profile of the atmosphere due to the 
changes in gas opacity. However, the overall spectra of mod- 
els will look quite different in regions where CO absorbs 
strongly, such as the 4.5 /im feature, so future, more com- 
prehensive fits of sulphide cloud models to observations will 
have to include non-equilibrium models. 

5. DISCUSSION 

5.1. Formation of Clouds 

Clouds must form in brown dwarf atmospheres as they cool; 
there is no way to avoid the condensation of different species 
as the atmosphere reaches lower effective temperatures. In 
these models, we parameterize the opacity of clouds by creat- 
ing a distribution of cloud material in the atmosphere which 
has a distribution of cloud particle sizes. Within the models, 
we can change those distributions. A cloud that sediments 
into a small number of large particles will settle into a thin 
layer and will not significantly change the emergent spectrum; 
the same cloud material organized into an extended cloud with 
small particle sizes will have a dramatic effect on the model 
spectrum. 

For these reasons, we require a model of cloud particle sizes 
and distribution as well as the underlying chemistry. When we 
consider models that include new or different clouds, we do 
not change any of the underlying chemistry of condensation; 
we change the opacity of the condensate particles and in doing 
so modify the effect that the cloud formation has on emergent 
flux. 

5.2. Sulfide or Silicate Clouds? 

Burgasser et al] (120101) invoked the reemergence of silicate 
clouds to explain the spectrum of Ross 458C. We suggest in- 
stead that the initial emergence of sulfide clouds would have 
a similar effect on the spectrum and provide a more natural 
explanation for the results. 

From observations, it is clear that the range in T dwarf col- 
ors just following the L/T transition is small; spectra of T 
dwarfs show no evidence that clouds still affect the emergent 
flux for objects slightly cooler than this transition. This ob- 
servation suggests that the iron and silicate clouds have dis- 
sipated between 1400 and 1200 K (for typical field dwarfs) 
and are no longer important in T dwarf atmospheres. If iron 
and silicate clouds were sometimes important in T dwarf at- 
mospheres, we would expect to see a population of relatively 
quite red objects at effective temperatures between that of 
Ross 458C and the L dwarfs; no brown dwarf with these prop- 
erties has been observed. 

As T dwarfs cool, the range in observed infrared colors in- 
creases; a population of red T dwarfs develops, which are red- 
der in the near-infrared than cloudless models predict. Based 
on these observations, we favor a mechanism that cannot 
strongly alter r e ff~900-1200 K T dwarf atmospheres, but nat- 
urally reddens r e ff<800 K T dwarfs. 

The emergence of sulfide clouds provides that natural ex- 
planation for this range in T dwarf colors at low effective tem- 
peratures. Just as the iron and silicate clouds condense in the 
photospheres of L dwarfs and change their observed spectra, 
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FIG. 12. — UGP S 0722-05 ne ar-in frared spectrum comparison. Two different models are compared to the observed spectrum of UGPS 0722-05 from 
ILucas et"al] {2010). As in Figure [TTI the left panels show the near-infrared spectra; the right panels show the same spectra and models with near- and mid- 
infrared photometry. Yellow points show J, H, and K photometry; orange show Spitzer [3.6] and [4.5] photometry; red show WISE Wi, W2, and W3 photometry. 
The filters for the photometric bandpasses are shown with corresponding colors along the bottom. The top plot shows a cloudless model spectrum that best 
matches the data, which has an effective temperature of 550 K and surface gravity log g=5.0. The bottom plot shows the best matching cloudy spectrum using 
sulfide clouds; it has an effective temperature of 500 K, log g=4.5, and /s C d=5. 



the sulfide clouds begin to condense in the photospheres of T 
dwarfs with temperatures cooler than 900 K, changing their 
observed spectra. 

We have not yet investigated whether the sulfide clouds will 
have identifiable spectral features that would confirm their 
presence in T dwarf atmospheres, but given the features in 
the sulfide indices of refraction (see Figure |2) these features 
would likely be in the mid-infrared. 

5.3. Outstanding Issues In T Dwarf Models 

There are several challenges in modeling T dwarfs that have 
not yet been addressed in these calculations. Because of the 
high densities in brown dwarf atmospheres, sodium and potas- 
sium bands at optical wavel engths are extremely pressure- 
broadened in T dwarf spectra (|Tsuii et al.ll999"l iBurrows et alj 
I200a lAllardet all 120051. |2007|). The wings of these broad- 
ened bands extend into the near-infrared in Y and J bands, 
creating additional opacity at those wavelengths. For these 
calculations, we use the line broadening treatment outlined in 



IBurrows et al.l (|2000j), which is somewhat ad hoc and poten- 
tially creates some inaccuracies in the model flux in Y and J 
bands. A calculation of the molecular potentials for potassium 
and so dium in these high pressu re environments, as is carried 
out in lAllard et al.1 (|2005, 2007), would improve the accuracy 
of these models. 

Another challenge in modeling T dwarf spectra is the inad- 
equacies of methane opacity calculations; methane is the only 
important gas-phase absorber with inadequate opacity mea- 
surements or calculations. Uncertainties in methane absorp- 
tion bands create inaccuracies in T dwarf mod els, especially 
in // band where it is a very strong absorber (ISaumon et al.l 
I20T1) . 

5.4. Breakup ofNa2S Cloud 

Sulfide clouds could form partial cloud layers with patchy 
clouds. One way to infer patchy cloud cover is to ob- 
serve variability in photometric colors; variability can indi- 
cate high-contrast cloud features rotating in and out of view. 
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Radigan et al. (2012) studied objects at the L/T transition and 
inferred that the iron and silicate clouds could be in the pro- 
cess of breaking up and forming patchy clouds in those atmo- 
spheres. A similar study of the variability of cool T dwarfs 
could reveal a similar physical process in sulfide clouds. 

5.5. Constraining Cloud Models with More Data 

A larger number of high fidelity spectra of the coldest T 
dwarfs would help to constrain these cloudy models. Cur- 
rently there are a few objects with effective temperatures 
cooler than 700 K with moderate resolution spectra. A larger 
sample of objects would give us better statistics on the over- 
all population of T dwarfs, with different surface gravities, 
metallicities, and cloud structures. 

5.6. Water Clouds 

At cooler effective temperatures, water clouds will con - 
dense in brown dwarf atmospheres dBurrows et al.l 12003b . 
Oxygen is 300 times more abundant than sodium in a solar- 
composition gas and the silicate clouds only use 20% of the 
total oxygen in the atmosphere, so water clouds will be much 
more massive and important in shaping th e emergent flux. 
As m i ssions like WISE fin d colder objects (Kirkpatric k et al.l 
1201 U Cushing et al.l 1201 ll) and these objects are observed 
spectroscopically, future models of brown dwarfs will need 
to include the condensation of these more volatil e clouds. 

Before the wa t er clo uds condense, lLoddersI (jl999); 
lLodders & Feglevl d2006l) predict that RbCl and CsCl will 
cond ense; however , the abundances of Cs and Rb are very 
low ([Lodders 2003) so these clouds will be optically thin. If 
equilib rium conditions prevail, NH4H7PO4 would al so con- 
dense ( Fegley & Lodders 1994; Vi sscher et al.l 120061) with a 
mass similar to that of the Na2S cloud. Whether NH4H2PO4 
condenses or P remains in the gas phase as PH3 (as on Jupiter 
and Saturn) deserves further study to examine potential effects 
on the spectra of the coolest brown dwarfs. 

5.7. Application to Exoplanet Atmospheres 

Observations and models of T dwarfs provide a testbed to 
study planetary atmospheres. While brown dwarfs are more 
massive than planets, the atmospheres of T dwarfs have sim- 
ilar effective temper a tures to those of young giant planets 
dBurrows et al l 119971: iFortnev et all 12008b . The study of T 
dwarfs provides crucial tests of cloudy atmosphere models 
that will be applicable to directly-imaged exoplanet atmo- 
spheres. 

Cloud models designed originally for brown dwarfs are al- 
ready being applied to exoplanets. Cloud models with iron 
and silicate clouds were originally developed to model L 
dwarf atmospheres; these models have been successfully ap- 
plied to observations of the only directly imaged multiple 
planet system, HR 8799. Several studies of the HR 8799 
planets have shown that the iron and silicate clouds play 
a significant role in their atmospheres (Marois et al.l l2008t 



B arm an et al 11201 It iGalicher et al.ll201 U lBowler et al.ll2010t; 
Currie et al l 1201 It iMadhusudhan et al.l 1201 ll iMarlev et all 
20121) . 



As instruments like the Gemini Planet Imager and SPHERE 
begin to discover new planets in the next few years, we will be 
able to apply brown dwarf models to colder planetary atmo- 
spheres in which clouds will likely play a key role in shaping 
their spectra. 

6. SUMMARY 

Cloud formation is a natural and unavoidable phenomenon 
in cool substellar atmospheres. At temperatures cooler than 
those of L dwarfs, chemistry dictates that additional conden- 
sates, beyond the silicates and iron, must form. We have ex- 
amined the effect of including the most abundant of these 
lower-temperature condensates, Cr, MnS, Na2S, ZnS, and 
KC1, in brow n dwarf model atmospheres. Within the frame- 
work of the Ackerman & Marie v (2001) cloud model, we 
have investigated the opacity of these clouds over a wide 
range of parameter space, across the relevant range of T 
dwarfs, to the warmest Y dwarfs. From our suite of mod- 
els from 400 to 1300 K, log g=4.0 to 5.5, f sed =2 to 5, we have 
shown the likely role that these low-r e ff clouds, dominated by 
sulfides, play in these cool atmospheres. 

Model spectra were compared to two T dwarfs, Ross 458C 
and UGPS 0722-05. These two objects have red near-infrared 
colors and are not well-matched by cloudless models. Model 
spectra that include the sulfide clouds match the observed 
spectra of both objects more accurately than cloudless mod- 
els. 

The photometric colors of the cloudy models were calcu- 
lated and compared to the full population of brown dwarfs 
with parallax measurements. This analysis shows that the sul- 
fide clouds provide a mechanism to match the near-infrared 
colors of observed brown dwarfs. The agreement is partic- 
ularly good in J-H, while in J-K the models are some- 
what too red. WISE observations of the coolest T dwarfs and 
warmest Y dwarfs indicate the these new models fit observa- 
tions better than cloud-free models. 

Our results indicate that understanding the opacity of con- 
densates in brown dwarf atmospheres of all T e ff is necessary 
to accurately determine the physical characteristics of the ob- 
served objects. 
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